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Abstract 

Some features of integrable lattice models are reviewed for the case of the six-vertex model. 
By the Bethe ansatz method we derive the free energy of the six-vertex model. Then, from the 
expression of the free energy we show analytically the critical singularity near the phase transition 
in the anti-ferroelectric regime, where the essential singularity similar to the Kosterlitz-Thouless 
transition appears. We discuss the connection of the six-vertex model to the conformal field theory 
with c = 1. We also introduce various exactly solvable models defined on two-dimensional lattices 
such as the chiral Potts model and the IRF models. We show that the six-vertex model has rich 
mathematical structures such as the quantum groups and the braid group. 

The graphical approach is emphasized in this review. We explain the meaning of the Yang- 
Baxter equation by its diagram. Furthermore, we can understand the defining relation of the 
algebraic Bethe ansatz by the graphical representation. We can thus easily translate formulas of 
the algebraic Bethe ansatz into those of the statistical models. As an illustration, we show explicitly 
how we can derive Baxter's expressions from those of the algebraic Bethe ansatz. 

Keywords: exactly solvable models, the Yang-Baxter equation, statistical mechanics, phase tran- 
sitions, the Bethe ansatz, the six-vertex model, integrable lattice models, conformal field theory, 
two-dimensional lattice 

1 Introduction 

We introduce the six-vertex model defined on a two-dimensional square lattice. We describe the 
model in detail, since it gives an important prototype of many solvable lattice models defined on two- 
dimensional lattices [1]. The transfer matrix of the six- vertex model generalizes the XXZ quantum spin 
chain which plays a central role among integrable quantum spin chains [2, 3]. The eight-vertex model, 
which generalizes the six-vertex model directly, may be considered as the most important exactly 
I solvable model in statistical mechanics [4]. Moreover, many mathematical theories such as the algebraic 
Bethe ansatz [7] and quantum groups [5, 6] are closely related to the six-vertex model. Starting from 
the six-vertex model, one may have a wide viewpoint on various physical and mathematical topics 
related to solvable models. There are quite a large number of topics related to exactly solvable models 
in physics and mathematics [1, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 
28, 29]. 

*In "Classical and Quantum Integrable Systems: Theory and Applications" (ISBN 07503 09598), Institute of Physics 
Publishing, (2003) http://bookmarkphysics.iop.org/ (Chap. 5 , pp. 113-151.) 
^e-mail deguchi@phys.ocha.ac.jp 
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We explain in §2 some features of the six-vertex model defined on a square lattice. We introduce 
the Boltzmann weights and the transfer matrix for the six-vertex model. We review a method for 
diagonalizing the transfer matrix, which is called the coordinate Bethe ansatz, and we show the 
expressions of the free energy per site in the ferroelectric, the anti-ferroclcctric, and the disordered 
phases, respectively [30, 31]. The disordered phase is gapless, while the ferroelectric and the anti- 
ferroelectric phases have gaps [32, 33]. We derive a critical singularity appearing at the phase transition 
from the anti-ferroelectric to the disordered phases. We review the calculation of the singular part 
of the free energy through the analytic continuation, as shown in Ref. [1]. The critical singularity 
is very weak and has the essential singularity similar to the Kosterlitz-Thouless transition. We have 
thus derived the KT-like singularity through exact calculation. After reviewing the finite-size analysis 
of conformal invariance [34, 35, 36, 37], we discuss that the massless phase of the six- vertex model 
is related to the conformal field theory of c = 1 which has U{1) symmetry. The c = 1 CFT has 
the critical line where critical exponents change continuously with respect to some parameter of the 
model [38, 39, 40]. There are quite a few papers on the finite-size corrections of the integrable models 
[41, 42, 43, 44, 45]. (For a review, see [46, 16, 47].) The critical line is also characteristic to the 
Tomonaga-Luttinger liquid [48]. The existence of a critical line was first discovered by R.J. Baxter 
through the exact solution of the eight- vertex model [49] . 

In §3, we review various integrable models in statistical mechanics [1]. We briefly introduce the 
Ising model [8, 50, 51, 52], the Potts models [53, 54, 55, 56], and the chiral Potts model [57, 58, 59, 
60, 61, 62, 63, 64, 65, 66, 67, 68] and then the eight-vertex model [49, 4, 69, 70] and IRF models 
[71, 72, 73, 74, 75, 76, 77, 78]. In §4, we solve explicitly the Yang-Baxter equations for the six- vertex 
model. We introduce the algebraic Bethe ansatz [7, 79, 80, 81, 16]. Here we show that the Yang- 
Baxter equation of the algebraic Bethe ansatz can be expressed by graphs. In §5, we discuss some 
mathematical theories associated with integrable models such as the braid group [82, 83] and the 
quantum groups [84]. 

There have been novel developments in the mathematical physics associated with the six-vertex 
model [84, 20]. The integrable vertex models associated with various Lie algebras have been obtained, 
which generalize the six- vertex model [85, 86]. The crystal basis of the quantum groups is derived 
from the mathematical analysis of the corner transfer matrix which is fundamental for calculating the 
one-point functions of the vertex models and IRF models [87, 88]. Through the vertex operators, 
correlation functions of the XXZ spin chain or the six vertex model are obtained [20] . The dynamical 
Yang-Baxter equation [89] and the elliptic quantum groups [90, 91, 92] have also been extensively 
discussed. In fact, we can derive the R matrix of the eight-vertex model systematically from the 
elliptic quantum group through the twists [92]. Furthermore, the correlation functions of the XXZ 
model calculated with the g-vertex operators have been rederived for large but finite chains through 
the algebraic Bethe ansatz with the Drinfeld twists [93, 94]. Here we note that the (7- vertex operator 
can be defined only on the infinite chain, while the algebraic Bethe ansatz with the Drinfeld twists 
can be applied to any finite chain. By taking the thermodynamic limit, it has been shown that 
the two approaches indeed give the same results. These papers indeed illustrate nontrivial physical 
applications of the Drinfeld twits. It is recently found that the symmetry of the six-vertex model is 
enhanced at some particular coupling constants: the transfer matrix commutes with the generators of 
the SI2 loop algebra for the six- vertex model at roots of unity [95]. 

Let us discuss some physical motivations for the six-vertex model. The exact solution of the six- 
vertex model was originally introduced for studying the statistical mechanics of ferro-electrics such 
as the residual entropy and the ferroelectric transitions [30, 31]. However, it seems that the physical 
motivation of the six-vertex model for ferro-electricity has decreased. On the other hand, there are 
many different physical applications of of the six-vertex model. Here we consider a few examples, 
the domain wall theory [96, 97], crystal growth [98, 99, 100, 101, 102, 103] and the thermodynamics 
of the XXZ spin chain through the quantum transfer matrix [104, 105, 106, 107, 108, 109]. The 
crystal growth on surfaces has been discussed by applying exact solutions of the six-vertex model 
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[98, 99, 100, 101, 102] and some extensions [103]. The free energy of the six vertex model gives 
the equihbrium crystal shapes [98, 99]. The finite-temperature thermodynamics of the XXZ spin 
chain has been studied extensively through the quantum transfer matrix, which is a version of the 
inhomogeneous six- vertex transfer matrix [105, 106, 107, 108, 109]. The quantum transfer matrix 
is obtained by regarding one direction of the square lattice as the imaginary time or the inverse 
temperature [104]. There have been considerable efforts to evaluate thermal quantities analytically or 
numerically. Several functional equations on the eigenvalues of the transfer matrix have been devised 
[108, 109]. Finally, we remark that a universal relation between the dispersion curve and the ground- 
state correlation length in quantum spin chains are discussed by using the exact solutions of the vertex 
models [110]. 

In §2, we employ mainly the notation of Baxter's textbook [1] except for the transfer matrix. In 
§4, however, we briefly show how the notation in statistical mechanics is related to the notation of the 

algebraic Bcthe ansatz or the Quantum Inverse Scattering Method. The graphical illustration should 
be useful. Finally in §5 we discuss many connections of the six- vertex model to several mathematical 
developments such as the quantum groups. 

2 Solvable vertex models 

2.1 The six- vertex model 
2.1.1 Ice rule 

Let us consider a square lattice as a model of two-dimensional ferroelectric crystal. Molecules are 
placed on the vertices of the lattice. Arrows are placed on the edges of the lattice, which correspond 
to directions of dipole moments of hydrogen bonds. As a crystal of hydrogen bonding, we may consider 
ice, i.e., the crystal of water molecules. In this review, however, we simplify the molecular background 
of the model (For instance, see [9]). We assume that the dipole moments defined on edges take only 
two values ±1. 

At a vertex in the lattice, there are four edges. There are possible 16 configurations of the four 
edges around the vertex since each of the edges takes two values ±1. 
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Figure 2.1: Configuration of polarizations around a vertex: a, /3, 7,(5. The Boltzmann weight is 
expressed by i(;(q;,/3|7, S). Positive directions of dipole moments are given by upward or rightward. 

Let symbols a, /3, 7, b denote the values of dipole moments around the vertex. Due to charge 
neutrality, they should satisfy the following condition 

a + ^ = 7 + 5 (2.1) 
For an illustration, let us consider the case when a;,/3, 7 and b are given by -|-1. Then, a and /3 give 
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+ 2 to the vertex, while 7 and d deprive + 2 from it, so that the net charge around the vertex is kept 
neutral: a + /3 — 7 — 5 = 0. 

There are only 6 configurations satisfying the condition. The other configurations that do not 
satisfy the condition are not allowed in thermal equilibrium. Thus, wc call the model the six-vertex 
model. The condition (2.1) is sometimes called ice rule, since ice as a crystal consists of water molecules 
connected by hydrogen bonding. 

We denote by 1 and 2 the values of polarization 1 and —1, respectively. The symbols 1 and 2 are 
useful for matrix notation. Let p denote the notation of ±1 and k 1 and 2. Then, they are related by 
the relation: k = 1 + {1 — p)/2. 



(1) 



(3) 



(5) 



(2) 



(4) 



(6) 



Figure 2.2: Vertex configurations satisfying the ice rule. They have the Boltzmann weights w{a, /3\^, 6) 
as follows: (1) w{l,l\l,l); (2) tf(2,2|2,2); (3) 2|2, 1); (4) w(2, 1|1,2); (5) u'(l,2|l,2); (6) 
■u;(2, 1|2, 1). The configurations (1) and (2) are for the weight a, (3) and (4) for b, and (5) and 
(6) for c. 



2.1.2 Boltzmann weights 

It is a key idea in exactly solvable models that we define the model by the Boltzmann weights not by 
the energies of configurations. Let us introduce the Boltzmann weights for configurations around a 
vertex. For a vertex configuration a,P,j,S, we denote by e{a,l3\j,S) the energy at the vertex. Then, 
the Boltzmann weight for a temperature T is given by 

w{a, 6) = exp(-e(a, 5)/kBT) (2.2) 

Under the ice rule, there are only six configurations allowed round a vertex. Here, it is assumed that 
the energy of a configuration violating the ice rule should be infinite. We denote by ej the energy of 
the jth vertex configuration shown in Fig. 2.2. 

Under no external field, the Boltzmann weights must be invariant when reversing all the polariza- 
tions simultaneously. Thus, we have e\ = 62, €3 = £4 and £5 = ee, when there is no external field. We 
denote the Boltzmann weights as follows. 

w{l,l\l,l) = w{2,2\2,2) = wi = a 
u;(l,2|2,l) = k;(2,1|1,2) =k;2 = 6 

u;(l,2|l,2) = w{2,\\2,l) = = c (2.3) 

The Boltzmann weights of the zero-field six-vertex model have essentially only two parameters. 
For instance, we may choose a/c and h/c. Note that the probability for the vertex configuration of a 
is given hj a/(a + h + c), which does not change by replacing a, h and c with pa, ph and pc. 
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Let us consider 7r/2 rotation of the square lattice. If we rotate vertex configuration (1) of Fig. 2.2 
by the angle 7r/2 in the counterclockwise direction, then, it becomes vertex configuration (4). Under 
the 7r/2 rotation, the weight a is exchanged with the weight b, while the weight c does not change. 

2.2 Pctrtition function and the transfer matrix 

Let us discuss the partition function of the system. We now set the boundary conditions. Here, we 
consider the periodic boundary conditions for the two-dimensional lattice. We take a product of the 
Boltzmann weights over all the vertices of the lattice, and sum up the product over all the allowed 
configurations of arrows on the lattice. 

^= E n <<^J'b,\cj,d,) (2.4) 

con fig j '.vertices 

The partition function of the square lattice can be formulated as the trace of the products of the 
transfer matrices. Let us define the transfer matrix r of the six-vertex model. The matrix elements 
of the transfer matrix r acting on N lattice sites are given by 

'^bi,'.'.'.',bN = E '^i<^i^ bi\ai, C2)w{c2, b2\a2, cs) ■■ ■w{cN,bN\aN,ci) (2.5) 

C1,...,CJV 



ai 
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Figure 2.3: Matrix elements of the transfer matrix T^^'"'f^ 



Under the periodic boundary conditions, the partition function Z^^^i of N x N lattice is given by 
the trace of the iV'th power of the transfer matrices: 



N 



,\ ai,...,aM I I 



E 

ai,...,aM ^ ^ oi,...,ajv 



N 
2N 



Here Aj denotes the eigenvalue of the transfer matrix r. 
The free energy per site / is given by 

f = -kBT\ogZ^^,/{NN'). 



(2.6) 



(2.7) 



In the thermodynamic limit: N,N — > oo, the free energy per site is given by the largest eigenvalue 
^max of the transfer matrix r. 



2.3 Diagonalization of the transfer matrix 

2.3.1 The Yang-Baxter relations for six- vertex model 



Let us consider three sets of the Boltzmann weights: {wi,W2, w^) = (a, b, c), (a', b' ,c), and (a , b" , c ). 
We denote by r and r the transfer matrices constructed from the sets of the Boltzmann weights 
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{a ,b' ,c) and (a" ,b" ,c"), respectively. If the three sets of the Boltzmann weights satisfy the Yang- 
Baxter equation 

^ w{a, 7|ai, 02 )«;'(/?, bs]-/, a3)w" {bi,b2\a, (3) 
= ^ w"{P,a\a2,a3)w'{bi,-/\ai,P)w{b2,b3\'y,a), (2.8) 

then the transfer matrices r' and r" commute. The derivation of the commutation relation is given in 
Appendix A. We note that a graphical presentation of the Yang-Baxter equation (2.8) will be shown 
in Fig. 4.1. 

Let us define the parameter A as follows 

. + 62 - c2 



2ab 



(2.9) 



For the zero-field six-vertex model, we can show that if the two sets of the Boltzmann weights have 
the same value of the parameter A, then their transfer matrices commute. We shall explicitly discuss 
in §4 that it is indeed derived from the Yang-Baxter equations (2.8). 



2.3.2 The coordinate Bethe ansatz 

Let us consider the matrix clement t^^'...^n^ of the transfer matrix r. Due to the ice rule, we may 
express the suffix ai, ... ,aj\i hy the positions of the value 2, as follows. Suppose that there are n suffices 
given by the value 2 among the N suffices ai , . . . , oat. The n suffices are expressed as , , ■ ■ ■ , o-xn 
where Xj's are in increasing order: xi < X2 < ■ ■ ■ < Xn ■ Then, the entry oi, . . . , oa? is equivalent to 
the set of the Xj's: xi, . . . ,Xn. For an illustration, let us consider the case N = 5 and n = 3. Then, 
{xi,X2,xs) = (1,3,4) corresponds to (01,02,03,04,05) = (2,1,2,2,1). Thus, the matrix element 
'^bi .ThM ^^"^ denoted briefly as Tyl^;;;^y^ . 

Let us now discuss how to solve the secular equation: rg = Ag. Here, the transfer matrix r is a 
2^ X 2^ matrix, g is a 2^-dimensional eigenvector with eigenvalue A. In terms of matrix elements, 
the secular equation is written as 

E ''yt'-t^aiyi, ■■■,yn) = ^gixi, . . . , x„) (2.10) 

yi,-,yn 

Here, g{xi, . . . , Xn) denotes the matrix element of the vector g for the entry [xi, . . . , Xn). In the coor- 
dinate Bethe ansatz, we assume the following form for the matrix element of the possible eigenvector 
9 

g{xi, ...,Xn)= XI i^Pi^i ^ kpnXn) (2.11) 

PeSn 

Here, Sn denotes the symmetric group of order n, and P is a permutation of n letters, 1,2, ... ,n, 
where P maps j into Pj. The expression (2.11) is called the Bethe ansatz wave-function. If the vector 
g whose elements are of the form (2.11) is an eigenvector of the transfer matrix, then we call it a Bethe 
ansatz eigenvector. 

For general n, the vector (2.11) becomes an eigenvector of the transfer matrix, if the wave-numbers 
kj's satisfy the Bethe ansatz equations. They are given by the following 

n 

exp{iNkj) = {-iy'-^l[exp{-ie{kj,ke)), for j = l,...,n, (2.12) 
e=i 

where @{p,q) is defined by 

1 - 2Af ''' + r 

exp(-ie(p, q)) = , \ ^ , (2.13) 



6 



For the solutions fej's to the Bethe ansatz equations, the eigenvalue A of the transfer matrix is given 

by 

A{ki, ...,kn) = a'' L{zi) ■ ■ ■ L{zn) + b"" M{z{) ■ ■ ■ M{zn) , (2.14) 
where Zj = exp[ikj) for j = 1, . . . ,n and the functions L{z) and M{z) are defined by 

^, , ab+(c^-b'^)z , a"^ - - abz 

a{a — bz) b{a — bz) 

When we discuss the spectrum of an integrable model through the coordinate Bcthe ansatz, we 
often assume that all the eigenvectors of the transfer matrix are characterized by the Bethe ansatz 
wave-function (2.11). However, it is not certain whether the assumption is valid or not. Thus, we 
have to check it by other methods. In fact, there arc several mimerical studies on the validity of the 
completeness of the Bcthe ansatz for some integrable models. 

On the other hand, there is no doubt on the mathematical structure of the Bethe ansatz wave- 
function. We can derive the expression (2.11) by the algebraic Bethe ansatz through the 'two-site' 
model [79, 16]. (There is an instructive notice in [HI].) It was shown that the matrix elements of 
the product of B operators acting on the vacuum are given by the Bethe ansatz wave-function (2.11) 
with /cj's being generic. 

The Yang-Baxter relation leads to not only the integrability of the six-vertex model but also the 
systematic construction of the eigenvectors. In fact, we shall see in §4 that the algebraic Bethe ansatz 
is solely based on the Yang-Baxter equation. 

2.3.3 An example of the eigenvector 

For an illustration, let us consider the eigenvector g for the case of n = 1: g{x) = Ae-K.-p{ikx). Through 
a direct calculation, we have 

N x-1 N 

Y.^y9{y) = 5^a--^-^6^-^+^-^c25(y)+ a^'+^-y-^by-^-'c^giy) 

y=l y=l y=x+l 

= (a^ L{z) + 6^ M{z)) g{x) + ^—^-^(1 - z"") (2.16) 
Therefore, g{x) becomes an eigenvector if the Bethe ansatz equation 2;^ = 1 is satisfied. 



2.4 The free energy of the six-vertex model 
2.4.1 Three phases of the six- vertex model 

There are three phases for the zero-field six-vertex model. They are given by the regions of the 
parameter A: ferroelectric phase when A > 1; anti- ferroelectric phase when A < — 1; disordered 
phase when — 1 < A < 1. It is found that the disordered phase (—1 < A < 1) is gapless (massless), 
while the ferroelectric phase (A > 1) and the anti-ferroelectric phase (A < — 1) are gapful (massive). 

Let us consider the phase diagram of the six vertex model shown in Fig. 2.4. We recall that the 
ratios a/c and 6/c determine the model. Here we note that the number of independent parameters 
is given by two, since the overall normalization factor is arbitrary. The regimes I and II give the 
ferroelectric phase, the regimes III and IV are anti-ferroelectric and disordered, respectively. In terms 
of the Boltzmann weights, the regime I is given hy a > b + the regime II by 6 > a -|- c, and the 
regime lY hy a + b < c . Then, the regime III is given bya<6-i-c, 6<c-|-a and c < a + b. 

Let us derive the three phases through an intuitive argument. For the ferroelectric regime of 
A > 1, we have a? + b'^ — c? > 2ab , which leads to the inequality: |a — 6| > c. When a > b, we have 
a > b + c. Thus, the configuration where all the Boltzmann weights are given by the weight a should 
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be the largest contribution to the partition function Z. In fact, when n = in eq. (2.14), we have 
A = + . When a > b, the free energy per site, /, is given by ei: / = — A:sTloga. When b > a, 
we have / = €3 , similarly. 

For the anti- ferroelectric regime of A < —1, we have the inequality a + b < c. Thus, the vertex 
configurations for c should be more favorable than those of a or b. In fact, it is shown that the transfer 
matrix has the largest eigenvalue when n = N/2. Furthermore, if we send |A| to infinity, all the vertex 
configurations should be given by those of c. The phase is thus called anti-ferroelectric. 

We may explain the reason why it is called anti-ferroelectric. Let us consider the configuration (5) 
of Fig. 2.2. Wc see that the arrow coming from the left goes upward, while the arrow coming from 
the right goes downward. If all the vertex configurations on the square lattice are given by (5), then 
the lines coming from South West to North East and the lines coming from North East to South West 
occupy the lattice alternatively. This gives the anti-ferroelectric order. 



b/c 




1 



Figure 2.4: Phase diagram of the six-vertex model: regimes I and II are ferroelectric, regime III is 
disordered, and regime IV is anti-ferroelectric. 



2.4.2 Parametrization of the Boltzmann weights 

It is nontrivial how to parametrize the Boltzmann weights. Recall that there arc two independent 
parameters for the six-vertex model. Thus, if we consider A as a parameter, there is only another one. 
As we shall see later, it is related to the spectral parameter. 

We recall that the phases of the zero-field six- vertex model are classified into the following: A < — 1, 
-1 < A < 1 and 1 < A. 
(1) Anti-ferroelectric phase 

For A < — 1, we define a real parameter A by 

A = -cosh A (0<A) (2.17) 
and we parametrize the Boltzmann weights as follows 

a = psinh( "^ 2 ^ '^ ' ^ ~ P ^™^( ^ ^ ^ ) ; c = psinhX, (— A < v < A) (2-18) 
where p is the normalization factor. Let us define the rapidity a for the wavenumber k as 

- e"*" sini(a - iX) ^ ^ 

eMik) = = ^ (2.19) 
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By replacing the wave numbers p and q with the rapidities a and /? in eq. (2.13), the phase factor 
0{p, q) is written as follows 



^2X_^-^ia-0) sini((a-/j)-2^A) 
exp(-^e(p,,)) = ^,,_.(._,) _ , = - ,,,i((,_^) + 2zA) ^'-'^^ 

(2) Disordered phase 

For — 1 < A < 1, we define a positive real parameter /x by 

A = -cos/x (0<Ai<7r) (2.21) 

and we parameterize the Boltzmann weights as follows 

. ,u — w, , , .u + w. . , , 

a = psm( — - — ), o = psm( — - — ), c = /9sm/x, {— /j, < w < /j,) . (2.22) 

Here p is the normalization factor. We define the rapidity a for the wavenumber k by 



sinh Ma — i/j.) 



eMik) = ^— = ^ (2.23) 

> e'ix+a_i sinhi(a + i/x) ^ ^ 

In terms of rapidities a and /3, the phase factor 6(p, g) is expressed as 

e^'^-e"-/^ sinhi((a-/3)-2i/x) 
exp(-ze(p,,)) = ^-^^^^ = - ,,,t,i((,_^) + 2^;.) ^'-''^ 

(3) Ferroelectric phase 

For A > 1, we define the real parameter A by 

A = cosh A (0<A) (2.25) 

and we may parameterize the Boltzmann weights as follows: When a > b (the regime I), 

a = psinh(^-^-^) , 6 = — psinh(^^-^-^) , c = psinh(A), {v < —X) . (2.26) 

When a < b (the regime II), 

a = -/9sinh( ^ ^ ^ ) , 6 = psinh C^ ^ ^ ) , c = /9sinh(A), {v > X) . (2.27) 

Here we recall that p is the normalization factor. The wavenumber k is related to the rapidity a by 

/,x e'^ — e~*" sini(Q;-iA) „„s 
^M^k) = , = — I7--TT (2-28) 



e 



;A-ia_i sini(Q; + a) 



The phase factor @{p, q) is written as 



g2A _ g-i(a-/3) sin 1 ( - /3) - 2z A) 



exp(-iG(p, g)) = ^rr— 7 — = 4 '- (2.29) 

V\ KF^HJJ g2A-»(a-/3) _ 1 sin i((a - /?) + 2zA) ^ ^ 

Let us consider curved lines given by changing the spectral parameter u; or u continuously. All the 
curves pass through the points (a/c, 6/c)=(l, 0) and (0,1). Except for the two points, however, any 
points on the horizontal axis: b/c = or the vertical axis: a/c = are never reached by the above 
parametrizations with finite values. 
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2.4.3 Expressions of the free energy 

(1) Anti-ferroelectric phase 

When A < — 1, the system is in the anti- ferroelectric phase. The free energy per site / is given by 

f = -keTloga- keT -—- + ^ -\ 2.30 

\ m=l / 

\ m=l / 

for —A < V < X. 

(2) Disordered phase 

When — 1 < A < 1, we have 

, r^-. , f°° sinh(u + ty)a: sinhfvr — u)a:: , 

/ = -kBTloga-ksT ^— 1 ^dx 2.32 

J-oo 2x smh TTX cosh /j^x 

, , , ^ f°° sinhfit — w)a:; sinhfTT — u)a; , 

= -kBTlogb-keT w—r^ r —dx 2.33 

J-oo 2x smh TTX cosh fj,x 

for — <w<ii. 

(3) Ferroelectric phase 

When a > b (the regime I), we have / = — AjfiTloga, while when b > a (the regime II), we have 
/ = -A:BTlog6. 

2.4.4 Correlation length 

Correlation length has been calculated for the six- vertex model [32, 1]. In fact, it is calculated for 
the eight-vertex model. In the ferroelectric and the anti-ferroelectric regimes of the six-vertex model, 
the correlation length is finite. It becomes very large near their phase boundaries to the disordered 
phase. In the disordered phase of the six- vertex model, the correlation length diverges throughout the 
regime. Thus, the disordered phase is critical. 

2.5 Critical singuleirity in the anti-ferro regime neeir the phase boundciry 

We now discuss a singular behavior of the free energy shown at the phase transition from the anti- 
ferroelectric to the disordered phases [1]. Here we note that the former one is gapful, while the latter 
is gapless. Recall that the phase boundary between the regimes III and IV is given by a + 6 = c. 
When T > Tg, the system should be in the disordered regime (— 1 < A < 1), while when T < Tc it is 
in the anti-ferroelectric regime (A < — 1). In the lower temperature, the system should be ordered. 

Let us calculate the analytic continuation of the high-temperature free energy (2.32) into the low- 
temperature phase, so that we can single out the singularity of the free energy near Tc and for T < Tc. 
First, we reformulate the integral of eq. (2.32) as follows. 

/■°° sinh(/L( -I- ?i;)a;sinh(7r — //)a:;^ ^ sinh(;Lt -I- ^ii)x cxp(7r — //)x^ (2 34) 

J-oo 2x sinh 7r.T cosh J-oo 2x sinh TTX cosh ^rc 

Here V denote the principal value integral. Second, wc set A to be a very small positive real number, 
and take a value v satisfying —A < i; < A. We consider the path in the complex //-plane: 

/x = Aexp(-z^), for O<0<7r/2. (2.35) 
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Along the path, we calculate the analytic continuation of the high-temperature free energy (2.32). 
Here, we also consider the path of w: w = vex.p{—i6) for < ^ < 7r/2. Then, we have 



^ roo — i sin(A + f )x exp(7r + iA)x^ 
2x sinh ttx cos Ax 



oo 



A + ^; ^ g-mA sinhm(A + i;) _ . ^ (-1)"^ cos ((m - l/2)7ri;/A) e-("-V2)^'A 
2 ^ mcoshmA ^ (m- 1/2) sinh ((m - l/2)7rVA) ^' ' 

The real part of the analytic continuation corresponds to the expression of the anti-ferroelectric free 
energy. Therefore, we obtain the singular part of the free energy 

_ ~ (-1)"- cos ((m - l/2)^t;/A) cxp (-(m - l/2)7r7A) 

fsing - il^Bi _ ^^2) sinh ((m - l/2)7r2/A) ^ '^'> 

We define the reduced temperature t by 

t = {a + h-c)/c. (2.38) 
In the low-temperature phase (T < T^) and near Tc, t is given by 

t~_l(A2_y2^ (2.39) 

8 

Approximately, t is given by i — A^/8. Near Tc, we have 

fsing ~ -4zfcBre-^'/^ cos(|^) (2.40) 

Thus, we have 

/ C( 

Near the free energy has an essential singularity. 

The singularity of the free energy is very close to that of the Kosterlitz-Thouless transition. In fact, 
calculating exactly, we can show that the correlation length ^ diverges at Tc as ^ oc exp (constant / ^/—t) , 
when T approaches Tc in the anti-ferroelectric phase. 

2.6 XXZ spin chain and the transfer matrix 

The logarithmic derivative of the transfer matrix of the six- vertex model gives the Hamiltonian of the 
XXZ spin chain. 

— logrL=_A = T~^—T oc Hxxz + constant (2.42) 
av dv 

where Hxxz is given by 

Hxxz = jj: {afaf^, + ajaj^, + Aafaf^,) . (2.43) 

3=1 

Intuitively, we may express it by tqv{v) ~ exp{—vHxxz)- 

The XXZ spin chain and the six-vertex transfer matrix have the same eigenvectors in common 
thanks to eq. (2.42). Taking logarithm of the Bethe ansatz equations, we have 

M 

Nkj = 27r/j -J2®ikj,ki), for j = 1, . . . , M . (2.44) 



/ constant \ , ^ , 

fsing oc exp ] (2.41) 
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where M is the number of down spins. (We assume 2M < N.) Here Ij is an integer if M is odd and 
half an integer if M is even. 

The ground state of the XXZ spin chain for A < 1 was obtained by Yang and Yang [3]. The 

ground state is specified by the integers Ij = j — (M + l)/2 for j = 1, . . . ,M. When N is very 
large, the distribution of fc^'s becomes continuous. The number of /cj's between k and k + dk can be 
approximated by N p{k)dk. Thus, we have the integral equation of p{k) 

27rp{k) = 1+1'^ ^^iMlp(k')dk' (2.45) 
J-Q ok 

where Q is determined by the normalization condition 

rQ 

/ p{k)dk = M/N . (2.46) 
J-Q 

The integral equation (2.45) can be solved by changing the variable k to rapidity a, and then by taking 
the Fourier transform for the half-filling case M/N = 1/2. When M/N is close to 1/2, the integral 
equation can be solved by the Wiener- Hopf method [3]. 



2.7 Low-lying excited spectrum of the transfer matrix and conformal field theory 

In this section, we assume that the low excited spectra of the transfer matrix of gapless models should 
be characterized by the conformal invariance, if the system size is large enough. The assumption is 
not rigorous, however, there are many studies which confirm it numerically for integrable models. We 
review the finite-size corrections for the XXZ spin chain and the six-vertex model. 



2.7.1 Finite-size corrections 

Let us consider a conformally invariant field theory defined in the two-dimensional Euclidean space 
with coordinates ri and r2- The energy momentum tensor Tf^^, for fj,,^ = 1,2 should be symmetric 
and traceless due to the conformal symmetry. Introducing the complex coordinates z = ri + ir2, 
z = ri — ir2, we define the chiral operator : T = (Tn ~ T22 — 2iTi2)/4, and the anti-chiral operator: 
T = (Til — T22 + 2iri2)/4. The operator T (or T) depends only on the variable z (or z). 
The energy momentum tensor has the operator product expansion 



T(,)Tfe) = ^^ + J^ + ^?^ + ... (2.47) 

[Zi - Z2r [Zl - Z2Y Zi - Z2 



Here c is called the central charge. We define the operators L„ by the expansion: T{z) = Z]^_oo " ^• 
The operator product expansion (2.47) corresponds to the Virasoro algebra 

[Ln, Lm] = {m- n)Lm+n + J^i^^ " rn)dm+n,o (2.48) 
Under a conformal transformation z ^ w, the energy momentum tensor is transformed as 

T{z)=(^^^'f{w) + ^{w,z} (2.49) 

where the symbol {-zi;, z} denotes the Schwarzian derivative: {d"'w / dz"^) / {dw / dz) - {'i/'i){d'^w/dz^Y /{dw/dz^ . 
Let us consider the conformal mapping from the z-plane to a cylinder of circumference L: 

2; — w = — log 2; (2.50) 
27r 
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Here w = T — ix with imaginary time t = it. The Hamiltonian H on the cyhnder is given by the space 
integral of the (1,1) component of the energy momentum tensor (Tcyi)^u 



H = dX {T^l {W) + fcyl {w)) 



= T^Lo + Lo)-^ (2.51) 
Here we have used (2.49). For the momentum operator on the cyhnder, we have 

27r 

= — (Lo-Lo) (2.52) 

Let us now discuss the apphcation of the formulas (2.51) and (2.52) to the quantum spin chains. 
We assume that the low excited energies should be gapless and conformally invariant. In other words, 

we assume that the excitations near the ground state have a linear dispersion relation. Let v denote 
the velocity of the linear dispersion. Then, for the ground-state energy Eq, we have 

TTVC 

Eo = Le^--^ (2.53) 
and for the excited energy E^x s-nd the momentum Pex 

we have 

'2tTTV — — 

Eex-Eo = —{h + h + N + N) 
2tt - 

Pex-Po = —{h-h + N-N) (2.54) 

where h and h are the conformal weights related to the zero modes of the field, and the eigenvalues of 
and N are given by non-negative integers. 

There is another viewpoint on the finite-size scaling. Let us consider the t axis as the space axis 
for an infinitely long quantum spin chain, and x axis as the imaginary time axis. Here we assume that 
L = v(3 = v/T. Thus, our system now becomes the quantum spin chain in the finite temperature T. 
Replacing £"0 with where / denotes the finite-temperature free energy of the spin chain, we have 
/ = —TrcT^/6v. Thus, we may calculate the specific heat C by the formula C = —Td'^f/dT'^, and we 
have 

C = gT (2.55) 

2.7.2 The free Boson: CFT with c = 1 

Let us consider a free Bose field ip{x, t) defined on a cylinder of circumference L. The Lagrangian is 
given by 

C- = \gj dx[{dt^f - {dx^f} (2.56) 

We define the mode ipn by the Fourier expansion ^p{x,t) = '/'n(^) exp(— 27rmx/L) . From the 
canonical quantization, we have the conjugate momentum 7r„ = gLif^n, and the commutation relation 
[<^n,7rm] = i^nm ■ With the operators a„ and a„ for n 7^ satisfying 

[an, am] = nSn+m [an, am] = [fln, Qm] = nSn+m , (2-57) 



the Fourier mode is expressed as (fn = i {an — a-n) / {n^/4:^^g) , for n / 0. The Hamiltonian is given 

by 

1 27r 

H = Yj/^ + 2Z ^ (o-nCln + a-nOn) (2.58) 
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Hereafter, we assume g = l/iir. The convention is consistent with the conformahy invariant partition 
functions. 

We now discuss the compactification of the Boson with radius R. Suppose that the field operator 
(f takes its value only on the circle of radius R. In other words, wc may identify (p with {p + 2itR. 
Then, the eigenvalue of the momentum ttq conjugate to ipo is given hy n/R for an integer n. Here we 
recall that the wavenumber of a one-dimensional system of size L is given by 27rn/L (n G Z), and also 
that the range of ipo is given by 2irR, which corresponds to L. Furthermore, we may assign on the 
operator ip the boundary condition for an integer m 

ip{x + L,t) = ip{x,t) + 2TrmR (2.59) 

Here, the integer m is called the winding number. The mode expansion of ip is given by 

cp{x, t) = ^o + —TTot + + ^ y - f a„e2'^^"(--*)/^ - a_„e2'^^"(-+*)/-^) (2.60) 

In terms of the coordinates z = exp(27r(r — ix)/L) and z = exp(27r(r + ix)/L), p{x, t) is given by the 
sum of the holomorphic and anti-holomorphic parts: (p{z, z) = (f){z) + ^{z) . Here, they are given by 

</'(^) = ?-i«olog(^) + iEr"fc^"' <^(^) = ^-molog(z) + i^iafcZ-^ (2.61) 

with ao = n/ R+mR/2 and oq = n/ R — mR/2 . Here we can show that the operator J{z) = id(p(z)/dz 
is the ^7(1) current operator. 

Making use of Noether's theorem, we have 

T{z) = -\ : {^f : , f{z) = -\ : {m)f ■ (2-62) 

Here : : denotes a proper normal ordering. Then, through the Laurent expansion of powers of z, we 
have 

-j^ oo 

^0 = 2^*0 + X! (^-n^n (2.63) 

rt=l 

Thus, the conformal weights hnm and hnm arc given by 
2.7.3 The XXZ spin chain and CFT with c = 1 

We discuss the finite-size corrections to the XXZ spin chain. The finite-size corrections to the ground 

state energy is calculated in [42, 43, 44, 45] based on the method with the Euler-MacLaurin formula 
[41]. (For a review, see [16, 46, 47].) The result is 

Ee. = Le^-^ + ^((Am^+(^ + iV + Ar) (2.65) 
27r 

Pe^-Po = 2kFAD + —{ADAM + N - N) (2.66) 

Here the term Coo denotes the ground-state energy per site. The Fermi velocity is obtained by the 
derivative of the dressed energy [16] with respect to the rapidity at the Fermi level. 

The central charge c is given by 1. AD and AM are integers. AM denotes the change in the number 
of down spins, and AD the number of particles jumping over the Fermi sea through the backscattering. 
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We note that the difference of the conformal weights (2.64) is given by hnm — hnm = nm. Thus, AM 
and AD correspond to n and m of the c = 1 CFT, respectively. N and N are derived from the 
particle-hole excitations near the Fermi surface. The Fermi wavenumber kp is given by kp = ttM/L 
where M is the number of down spins. If the dispersion is linear, kp is consistent with the number of 
particles M. 

The parameter ^ is given by the dressed charge, which is defined by an integral equation. We note 
that the sum of the conformal weights (2.64) is given by hn,m + ^n,m = ri^/R^ + m^B^/A. Thus, the 
dressed charge ^ corresponds to the radius R of the c = 1 CFT, ^ = R/2. Under zero magnetic field, 
the dressed charge ^ or the radius R is given by [33] 

/ 27r \V2 

R = (0 < M < tt). (2.67) 

Vtt - n) 

3 Various integrable models on two-dimensional lattices 

3.1 Ising model and Potts model 

3.1.1 Ising model 

Let us consider the Ising model defined on a square lattice. Each lattice site has spin variable which 
takes the two values ±1. We denote by Oj the spin variable of lattice site j. The Hamiltonian of the 
Ising model is given by 

n = -JY.aiaj (3.1) 

where the symbol denote that sites i and j are nearest neighbors, and we take the sum over all 
the pairs of adjacent sites on the lattice. 

There have been many papers written on the two-dimensional Ising model [8]. However, it should 
be remarked that correlation functions are calculated exactly for the two-dimensional Ising model in 
the scaling limit [50, 51]. 

We remark that exact solutions are discussed for the Ising model defined on various two-dimensional 
lattices such as the Kagome lattice (For a review, see Ref. [52] ). 

3.1.2 Self-dual Potts model 

The Potts model generalizes the Ising model into a p-state model with p > 2 [1, 53, 54]. Let us consider 
the three-state Potts model defined on a square lattice [53]. Each lattice site has spin variable which 
takes three values 1, 2, 3. We denote by the spin variable of lattice site j. The Hamiltonian of the 
Potts model is given by 

n = -JY.S{ai,aj) (3.2) 

Here the symbol (5(a, 6) denote the Kronecker delta. 

In general, the Potts model is not solvable. However, at its criticality, it is equivalent to some 
variant of the six- vertex model, and is solvable. Thus, the self-dual Potts model is integrable [55, 56]. 

3.1.3 Ashkin-Teller model 

With each site i we associate two spins: Sj and ti. They take values ±1. The Hamiltonian of the 
model is given by 

= - ^ \K2{siSj -\- titj) + KiSiSjtitj] (3.3) 

It is known that the Ashkin-Teller model and the eight-vertex model are in the same universality class 
[38] . The universality class is described by the c = 1 CFT with the twisted boson [40] . 
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3.2 Chiral Potts model 



3.2.1 General case 

There is another version of the Potts model [1, 53]. We may assign the chiral symmetry, or Z/pZ 
symmetry on the Potts model, which we call the chiral Potts model. 

We now explain the most general chiral Potts model defined on a square lattice [64]. Let a and 
b denote the spin variables defined on two nearest-neighboring sites. The interaction energy between 
the spins depends on the difference n = a — b (mod N) as 

N-l 

S{n) = EjU^"" u = e^^l^ (3.4) 

We note that £{N + n) = £{n). The parameter Ej can be written as 

E- A N 

^ = -K,cv^., for J = !,...[[-]] (3.5) 

where Kj and Aj constitute — 1 independent variables, and the symbol [[•]] denotes the Gaussian 
symbol. For a real number x, [[x]] denotes the biggest integer not larger than x. We also assume that 
En-j is complex conjugate to Ej-. E^^j = Ej. When N is odd, we have 

S(n) [[(^-1)/'" /27r \ ((-1)^ + 1) 

-|g= E 2i^, cos (|(jn + A,)) + K^M-ir \ ^ ^ (3-6) 

3.2.2 Integrable chiral Potts model 

Let us discuss the integrable restriction of the most general chiral Potts model defined on the square 
lattice. It has horizontal and vertical couplings. Suppose that spin variables a and b are located on 
two neighboring sites connected by a horizontal line. When the line goes rightward from a to b, the 
horizontal coupling has the energy £pq{a — b) and the Boltzmann weight Wpq{n). Here p and q are 
'rapidity' parameters. For spin variables c and d located on two neighboring sites connected by a 
vertical line, the vertical coupling has the energy £pq{c — d) and the Boltzmann weight Wpq(c — d) 
when the vertical line goes upward from c to d. 

The model is called solvable if the Boltzmann weights PF's and W's satisfy the star-triangle equa- 
tion 

N 

J2 Wqr{b - d)Wpr{a - d)Wpq{d - c) = RpqrWpq{a - b)Wpr){b - c)Wqr{a - c) (3.7) 

d=l 

The solution to the above equation is given by 

Wpqin) = ^P.(0) n L/^. • ^f^f;^) (3.8) 



with a constant Rpqr depending on the three rapidity variables. The constraint that the Boltzmann 
weight should have periodicity modulo N: Wpq{n + N) = Wpq{n) gives for all rapidity pairs p and q 

/ ,, \ ^ yN _ N N _ N 
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We can define k and k' such that 



k' 1 - fe< ..... 

x^ + y^ = k{l+x^y^) (3.11) 

where k"^ + (k'Y' = 1. Thus, the rapidities are placed on a curve of genus ^ > 1 of Fermat type. 

We make some comments. Multiplying the two equations of (3.9) and noting that p and q are 
independent, we can show eq. (3.10) and then eq. (3.11). The star-triangle equation is proven 
illustratively in the appendix of Ref. [60]. We note that the vertex- type formulation of the chiral 
Potts model is related to the tetrahedron equation [63] . The integrable chiral Potts model generalizes 
the self-dual Zjv model given by Fateev and Zamolodchikov [65]. An elliptic extension of the self-dual 
Zf^ model is introduced in Ref. [66], and the Yang-Baxter equation is proven for the model in Ref. 
[67]. Some nontrivial connections among higher-rank chiral Potts models, elliptic IRF models and 
Belavin's Zj<j symmetric model are explicitly discussed in Ref. [68]. 

3.3 The eight-vertex model 

Let us explain the eight-vertex model which was solved by Baxter in 1972 [4]. The Boltzmann weights 
w{a^ /3|7, (5; u) of are nonzero if the charge is conserved modulo 2: a + (3 = ^ + 5 (mod) 2. We have the 
six configurations around a vertex shown in Fig. 2.2 and the two in Fig. 3.1 . We assume that there 
is no external field. The weights therefore become symmetric: ei = €2, €3 = €4, €5 = ee, and ey = eg. 



(7) 



(8) 



Figure 3.1: Vertex configurations for the eight-vertex model. (7) 'u;(2, 2|1, 1); (8) u;(l,l|2, 2). 



iu(l,l|l,l) = w{2,2%2) = wi = asv 

u;(l,2|2,l) = w{2,l\l,2)=W2 = h^v 

u;(l,2|l,2) = w{2,l\2,l) = = c^v 

u;(l,l|2,2) = w{2,2\l,l) = Wi = d^v (3.12) 

We give a parametrization of the Boltzmann weights. We define the theta function 

00 

e{z-T) = 2pV4sin7r2 ]J(i _p2n-)^j _^2ng^p,^27rzz))(l -p2"exp(-27riz)) , (3.13) 

n=l 

where the nome p is related to the parameter r by p = exp(7rir) with Im r > . We also 
define the theta functions ^o(^) and 9i{z) satisfying 0a{z -|- 1) = {—l)°'9a{z) and 9a{z + t) = 
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ie '"^^^^'^^'^^Oi-aiz) for a = 0, 1 , where we define 6i{z) by 6i{z;t) = 9{z;2t). The Boltzmann 
weights asv{z),bsv{z),csv{z) and dsviz) are expressed as 

eo{z)eo{2r^) e,iz)eo{2r^) 
_ go(^)gi(2r?) ^1(^)^1(2^) 

- -^,(.-2,)^o(0)' '«^^^)--^o(.-2,)^o(0)- ^'-''^ 

3.4 IRF models 

3.4.1 Unrestricted 8V SOS model 

We introduce unrestricted Solid-on-Solid models. They are also called Interaction-Round-a-Face mod- 
els or IRF models, briefly [1]. 

To each site i of a two-dimensional square lattice, a spin a, is associated. Let i,j,k, and £ be 
the lattice sites surrounding a face (or a square), where i,j,k,£ are placed counterclockwise from 
the southwest corner. We assume that an elementary configuration is given by that of the four spin 
variables around the face, and the probability of having 0^,0^,0^,0^ is denoted by the Boltzmann 
weight w{ai,aj, Ofe, o^; z). Here the variable z is called the spectral parameter. 

For unrestricted 8V SOS model, Oj can take any integer. When sites i and j are nearest neighboring, 
then the states Oj and aj are said to be admissible if and only if \ai — aj| = 1 [69, 1]. The Yang-Baxter 
equations are given by 

^ w{a, b, g,f;z- w)w{f, g, d, e; z)w{g, b, c, d; w) 

9 

= ^w{f,a,g,e;w)w{a,b,c,g;z)w{g,c,d,e;z -w) , (3.15) 
9 

where the summation of the variable g is taken over all the admissible states. 
The Boltzmann weights are given by 

9(2n — z) 

w{d+l,d + 2,d+l,d;z,wo) = w{d,d - l,d - 2,d - 1; z) = ^ ' 



e{27j) 

w{d — l,d,d + l,d; z,wo) = w{d + l,d,d — l,d; z,wo) 



9{z) ^9{2r]{d + 1) + wo)9{2r}{d - 1) + wq) 



w{d + l,d,d + 1, d; z, wq) 



9{2r]) 9{2rid + wo) 

9{z + 2r]d + wo) 
9{2i]d + wo) 



1 ^ ^ 1 ^ ^ i9(z - 2rid - Wq) 

w{d- l,d,d- l,d;z,wo) = — (3.16) 

9{2rid + Wq) 



3.4.2 RSOS models 

Let us explain restricted Solid-on-Solid models (RSOS models) [71]. Let s denote the number of 
elements in S. Consider a s x s matrix C satisfying the following conditions [76, 77]: 

(i) Cab = Cba = or 1 

(ii) Caa = 

(iii) For each a £ S, there should exist b £ S such that Cab = 1 

For such choice of C, we impose a restriction that two states a and b can occupy the neighboring 
lattice sites if and only if Cab = 1- We call such a pair of the states (o, b) admissible. For the case of 
unrestricted models, the infinite matrix C satisfies the conditions (i), (ii), and (iii) with an infinite set 
S. 
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For an illustration, let us consider the restricted eight-vertex Solid-on-Solid model (the restricted 
8V SOS model), which we also call the ABF model [71]. For the A'^-state case, we have S = 
{1,2,..., N}. The nonzero matrix elements of C are given by Cjj-^-i = Cj-^-ij = 1 fov j = 1,2, . . . , N — 
1; other matrix elements such as Ci.iv and Cjv^i are given by zero. Setting wq = 0, then we have the 
Boltzmann weights of the ABF model. The Boltzmann weights satisfy the Yang-Baxter relations 
(3.15) with the finite set: S = {1,2,..., N}. 

Let us explain CSOS models, another type of RSOS models [75, 76, 77]. Here we assume that 
2Nr] = mi, where integer mi has no common divisor with A^. If wc set wq ^ 0, then we have the 
Boltzmann weights of the cyclic SOS model (CSOS model). Wc can show that the Boltzmann weights 
satisfy the Yang-Baxter relations with the finite set: S = {1,2, ...,A^} and the cyclic admissible 
conditions: Ci^jv = Cn,i = 1- 

We remark that the connection between the 6j symbols and the Boltzmann weights of IRF models 
is first discussed in Ref. [78]. 

3.4.3 Fusion IRF models and ABCD IRF models 

The 8VS0S model was generalized into the fusion IRF models. [72] IRF models associated with An ^ 
Lie algebra are constructed. [73] The IRF models associated with B^^^ C^^^ D^^^ typed Lie algebra 
arc also obtained [74]. 

In the IRF models, the one-point function, which is the magnetization per site, can be calculated 
by the corner transfer matrix method invented by Baxter [1]. 

3.4.4 Gauge transformations 

It is sometimes convenient to employ a gauge transformation 

w{a, b, c, d; z) w{a, b, c, d;z) — . (3-17) 

9a 

The transformed Boltzmann weights also satisfy the Yang-Baxter relations (3.15). For instance, we 
may set ga = exp(7rza/2) ^6{2r]a + wq) (a G Z). 

4 Yang-Baxter equation and the algebraic Bethe Ansatz 

4.1 Solutions to the Yang-Baxter equation 

4.1.1 Derivation of a solution for the six- vertex model 

Let us solve the Yang-Baxter equation for the six-vertex model. We recall that it is given by 

^ w{a, j\ai,a2)w' {P, 63I7, 03)11;" (61, 62IQ;, (3) 

= ^ w" {I3,a\a2,a3)w' {bi,j\ai,P)w{b2,b3\'y,a) . (4.1) 

The Yang-Baxter equation is illustrated in Fig. 4.1. 

There are 2^ x 2^ = 64 cases for the entries (oi, 02, 03) and (61, 62, ^'3)- Due to the ice rule, however, 
the Yang-Baxter equation is trivial unless 04 + 0,2 + 0,3 = 61 -|- 62 + ^'s- We have thus only 20 entries: 
(3C0)' + isCif + (3C2)' + (3C3)' = 1 + 9 + 9 + 1 = 20. 

Let us consider symmetries of the Boltzmann weights given in (2.3). If we exchange 1 and 2, the 
Boltzmann weights do not change. Thus, we have reduced 20 cases into 10 cases. 

(ai, a2, 03; 61, 62, ^3) = (1, 1, 1; 1, 1, 1), 
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(1,1, 2; 1,1,2), (1,1, 2; 1,2,1), (1,1,2; 2, 1,1), 
(1,2, 1; 1,1,2), (1,2,1; 1,2,1), (1, 2, 1;2, 1,1), 
(2, 1,1; 1,1, 2), (2, 1,1; 1,2,1), (2, 1,1; 2, 1,1) 

We now recall that the Boltzmann weights (2.3) have the symmetry 

w{a,l3\j,6) =w{j,d\a,l3) =w{l3,a\6,j) 



(4.2) 



(4.3) 



Combining these symmetries we can show that the Yang-Baxter equation for the two entries (1) and 
(2) are equivalent: (1) (ai, Q^s! ^i, ^2; ^3) ; (2) (63, 62, ^i; 03, 02; Qi)- Precisely, the L.H.S. (or R.H.S) 
of the Yang-Baxter equation of the case (1) corresponds to the R.H.S. (or L.H.S) of eq. (4.1). Thus, 
we have the three cases as follows. 



(1, 1,2; 1,1,2),(1, 1,2; 1,2,1), (1,2,1; 1,1,2) 
For the three cases, the Yang-Baxter equations are given by 



(4.4) 



ac a 



be b + ca c 

ab c = ba c + ccb 

I II I II , ' " 

cb a = ca b + be c 



(4.5) 



A nontrivial solution (a" , 6" , c" ) exists only if the determinant vanishes 

-he 



ae —oe —ea 
ec ba — ab' 
eb —ea —be 



We define the parameter A by 



abea b e 



{ar + {b' 



a^ + b^ 



a b 



ab 



A = 



+ 6^ - (? 
2d) 

The condition that the determinant vanishes is given by 

A = A' 



(4.6) 



(4.7) 



(4.8) 



Thus, the transfer matrices r and r commute, if the two sets of weights (a, 6, c) and {d ,1) ,c) have 
the same A. 
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In terms of the spectral parameter u, we can parametrize the three Boltzmann weights a, b and c. 
We express the weight as w{a, /3|7, 6; u). Let u and v be arbitrary. We denote w{a, /3|7, 5), w (a, /3|7, 6) 
, and w (a,/?|7, (5) as w{a, (3\^,S;u), w{a, (3\^,S;u + v), and w{a, Pl^yjS-jv), respectively. The Yang- 
Baxter equations are depicted in Fig. 4.2. As a solution, wc may have {a,b,c) = (psinh(u + 2rj),= 
psinhu, = psinh2?7 . Here we set A = cosh2r/ . The transfer matrices t{u) and t{v) commute: 
t{u)t{v) = t{v)t{u). 



4.1.2 Gauge transformations for vertex models 

Let us suppose that u;'s satisfy the Yang-Baxter equation. Then we can show that transformed weights 
u;'s defined by 

w{a, /3|7, 6; u) = (e)""'''^ cxp(K(a + 7 - /3 - 6)u)w{a, /3|7, 5; u) (4.9) 

also satisfy the Yang-Baxter equations [82, 83]. Here e = ±1, and the number n is arbitrary. 

The gauge transformation is important in the derivation of the Jones polynomial from the sym- 
metric Boltzmann weights of the six- vertex model under zero field [82]. It is also quite useful when 
we discuss the relation of the six-vertex model to the quantum group, as we shall see in §5 (see also 
the appendix of Ref. [95]). 



4.2 Algebraic Bethe ansatz 
4.2.1 R matrix and the L operator 

Let us diagonalize the transfer matrix by the method of the algebraic Bethe ansatz. [80, 81, 16] 

Let us introduce the notation of the matrix tensor product. We define the direct product A® B 
of matrices A and B. Let A^^^ denote the matrix element for the entry of column j and row k of the 
matrix A. Then, the matrix element of column (^1,^2) and (fci, ^2) is defined by 



(A< 



(4.10) 



We now define the R matrix of the XXZ spin chain. The element corresponds to the entry of 
column (a, b) and row (c, d). 



R{z) 



( R{zfA 
i?(z)lf 



R{z)\\ 

R{z)l\ 
R{z)ll 



R{z)\\ 

R{z)ll 
R{z)l\ 

Riz)ll 



Here a{z), b{z) and c{z) are given by 



Riz)ll \ 
Ri^W2 
R{z)l\ 
R{z)ll I 



( a{z) 



V 





c{z) 
b{z) 






b{z) 

c{z) 

a{z) ) 



(4.11) 



a{z) = sinh(z -|- 2ri) , b{z) = sinhz , c{z) = sinh2r7 . 



(4.12) 



Here, the functions a{z), b{z) and c{z) are equivalent to the Boltzmann weights a, b and c in §2 

We now introduce L operators. We write the matrix element the L operator with entry [j, k) as 



iLn{z)).k or Lni 



The L operator for the XXZ spin chain is given by 



sinh {zIn + Wn) 
sinh 277 cr+ 



sinh 2rj 
sinh {zIn - Wn 



(4.13) 



Here In and (n = 1, . . . , L) are acting on the n th vector space Vn- The L operator is an operator- 
valued matrix which acts on the auxiliary vector space Vq. The symbols denote (T+ = £^12 and 
a~ = E21, and a^,a''^, are the Pauli matrices 
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(1) 

b 



OLr. 



z-V a 




(2) 




Figure 4.2: (1) Ln{z)tU,f,^ ; (2) R{z)ll^^^^ 

In terms of the R matrix and L operators, the Yang-Baxter equation is expressed as 

R{z - t) {Ln{z) ® = {Ln{t) ® Ln{z)) R{z - t) (4.14) 

Here the tensor symbol in Ln{z) Ln{t) denotes the tensor product of auxihary spaces. 

The Yang-Baxter equation (4.14) gives the relation between the two products of 4 x 4 matrices. 
For an illustration, we consider the L.H.S. of (4.14). 

= ^i?(z-t)-;-(L„(.)0L„(t))^-^^ 



[R{z - t) ■ Lniz) Ln{t) 



61,62 



Cl,C2 



(4.15) 



C1,C2 



Here the symbol x denotes the product of matrices acting on the n th space Vn- Expressing the 



t — rj 




ai 



a2 



Figure 4.3: The Yang-Baxter equation: R(z - t){L{z) ® L{t)) = {L{t) ® L{z))R{z - t), eq. (4.14). 
operator products in the rath space Vn, the L.H.S. of (4.14) is written as follows 



[R{z - t) ■ Ln{z) (g) Ln{t) 



ai,a2 



61,62 



an,/?„ = E ^(^ - *)ci\'c2^ {Ln{z) ® urns \-n,Pn 



Cl ,C2 



= EE^(^-*)ci\'c2'^-(^)6iUn,7n^n(i)g|,„,/3„ (4.16) 



C1,C2 7n 



4.2.2 Monodromy matrix and the construction of the eigenvector 
We define the monodromy matrix by the product of the L operators (see also Fig. A.l) 

T{z) = Ln{z) ■ ■ ■ L2{z)Li{z) . 



(4.17) 
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The transfer matrix tqv{z) of the six vertex model is given by the trace of T(z) 

T,y{z)=tiT{z)=A{z) + D{z), where = ( ^[^| ) (4-18) 

The Yang-Baxter equation (4.14) leads to the commutation relation: R{z — t){T{z) (8) T{t)) = {T(t) ® 
T{z))R{z — t), from which wc have many relations among the operators A, B, C, and D. For instance, 
we have B{z)B{t) = B{t)B{z). Furthermore, we have 



We define the 'vacuum' by 



N 



|0) = |T)i|T)2---|T)iv (4.21) 
Multiplying A{z) and D{z) on the vacuum, we have 

A{z)\Q) =a{z-'nf\id) , D{z)\0) =h{z-r]f\Q) . (4.22) 

Let us consider the vector generated by the product of B operators. 

\M) = B{ti)---B{tMm (4.23) 

Then, through the commutation relations such as (4.19) and (4.20) wc can show [7, 16] that the vector 
\M) gives an eigenvector of the transfer matrix t%v{z) if rapidities ti,t2,----, tM satisfy the set of 
equations 

/ a{tj -ri) \^ _ c{tk - tj) bjtj - tk) ^ ( a{tj-tk) bjtk - tj) \ ^ 

\ b{tj -v)J b{tk - tj) citj - tk) . V bitj - tk) aitk -tj)) ' y . ) 

They are the Bethe ansatz equations (2.12) with different parametrization. For a set of solutions, 
ti,t2, ... ,tM io eqs. (4.24), the eigenvalue of the transfer matrix Tev{z) is given by 

= sinh^(. + ^) n ""^^t" " ^ + - ^) n " ' • 

/=i smh(ij-z) f^^ smh(tj-z) 

(4.25) 

4.2.3 Connection to the coordinate Bethe ansatz result 

Let us compare the result of the coordinate Bethe ansatz in §2. We consider the disordered phase: 

— 1 < A < 1. First, we change the variables w and a defined in §2 into u and ^ by it; = 2m — /x and 
a = 2C,— i/j,, respectively. Thus, from the expressions (2.22) we have (a, b, c) = (sin(/x — tt), sin u, sin /x) 
= ■j(sinh(— Z|Lt + zu), sinh(— iit), sinh(— i/Li)) and 

jl^Baxterr^A ^ ab + (c^ - b'^)zj ^ sinh((aj -iw- 2i^)/2) ^ sinh(-(Cj - iu)+ifi) 
a{a — bzj) sinh((Q;j — iw)/2) smh(Q — iu) 
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Here we recall that the symbol 1,^"-^^'^^ (^z) has been defined in eq. (2.15) in order to denote the 
eigenvalue of the transfer matrix of the six-vertex model. The expression (4.26) can be derived from 
the formula (4.25) of the algebraic Bethe ansatz as follows. First, we take the gauge transformation: 
b{z) —b{z), which corresponds to the case e = — 1 and k = in eq. (4.9). Then, we replace the 
variables z, ij's, and 2r] in (4.25) hy z — r] ^ —iu, tj — r] ^ —(j, and 2r] — > i/i, respectively. We have 
the following: 

a{tj — z) a{tj — z) smh.{tj — z + 2ri) sinh(— (Cj — m) + i/[i) (4 27) 

b(tj — z) b(tj — z) smh.{tj — z) sinh((^j — iu) 

and sinh^(2; + ??) ^ smh^ {—iu + ifi) . Thus, the formula (4.25) reproduces the expression (2.14) for 
the eigenvalues of the transfer matrix except for the normalization factor i^ . 



5 Mathematical structures of integrable lattice models 
5.1 Braid group 

5.1.1 The Yang-Baxter equation in operator formalism 

The Yang-Baxter equation in §2 gives a sufficient condition for the existence of commuting transfer 
matrices. However, there are other viewpoints on the Yang-Bax;ter equation. 
Let E-j^ denote the matrix given by 

{El)l = Sa,jS,,k a, 6 =1,2. (5.1) 

We define operators Xj(«)'s by 

(g>(j-l) ^(N-j-l) 

Xj{z)= ^ w{a,b\c,d;z)I®---®I®El®E^®I®---®I for j = l,...,N-l. (5.2) 

Here, the symbol w{a^h\c^d'^z) corresponds to and the essential part of Xj{z) is given by 

X{z) = ^ w{a, 6|c, d- z)El ® E^ (5.3) 

which is equivalent to R{z). In terms of the Xj(2;)'s, the Yang-Baxter equation is expressed by the 
following 

Xj{z)Xj+^{z + t)Xj{t)=Xj+^{t)Xj{z + t)Xj+^{z), for j = 1, . . . , iV - 1 . (5.4) 



5.1.2 The braid group 

The braid group for N strings is an infinite group which is generated by the generators 5i, . . . , 6Ar_i 
satisfying the defining relations 

bjbj+ibj = bj+ibjbj+i 

bibj = bjbi for |? — i| > 1 (5-5) 

Let us assume that the limit: \imz-^oo Xj{z) exists. Then, eq. (5.4) becomes 

Xj{oo)Xj+i{oo)Xj{oo) = Xj+i{oo)Xj{oo)Xj+i{oo) , for j = 1,...,N -1. (5.6) 
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This is nothing but the defining relations of the braid group. Thus, the Boltzmann weights of solvable 
models expressed in terms of the spectral parameter lead to representations of the braid group. 

We now show that from a given exactly solvable model, one can derive two different representations 
of the braid group [82]. Here, the gauge transformation (4.9) plyas a central role. This technical point 
is quite fundamental when we make connections of exactly solvable models to quantum groups and 
the Temper ley-Lieb algebra (for instance, see the appendix of Ref. [95]). 

We first consider the Boltzmann weights of the zero-field six- vertex model given by eqs. (4.12). 
Taking the infinite limit to them, we have the following 



lim X{z)/smh{z + 2ni) 



(I \ 

exp(-2r/) 

exp(-2r7) 

V 1 / 



(5.7) 



The representation of the braid group (5.7) leads to a link polynomial equivalent to the linking number. 

Let us now apply the gauge transformation (4.9) with e = 1 and «; = 1/2 to the weights given by 
eqs. (4.12) [82] . Then, from the transformed Boltzmann weights we have 



lim X{z)/ sinh(2; + 2ri) 



( 1 



cxp(— 2?;) 
V 



\ 

cxp(— 27/) 

1 - exp(-4?7) 

1 j 



(5.8) 



This matrix representation of the braid group leads to the Jones polynomial with q = exp(277) [82]. 

We remark that based on the representations of the braid group which are derived from the 
Boltzmann weights of exactly solvable models, we can construct various invariants of knots and links 
(see, for reviews, Refs. [83, 112, 113]). 



5.2 Quantum groups (Hopf algebras) 

From the quantum groups (the Hopf algebras) we can systematically construct representations of the 
braid group such as derived from the six-vertex model. Almost all the solutions of the Yang-Baxter 
equations can be constructed in some framework of quantum groups. Furthermore, the connection of 
solvable models to the quantum groups is useful for investigating nontrivial properties of integrable 
models. 

Let us introduce the quantum group Uq{sl2) , which is a g-analog of the universal enveloping 
algebra of 5^2. Generators X^,H satisfy 

KX^K-^ = q^^X^ , [X+, X-] = (5.9) 

q-q 

We may express K SiS K = q^. Taking the limit of q to unity, the relations are reduced into the 
commutation relations of sl2- 

The tensor product is defined by the following 

A{K) = K®K 
A(X+) = X+®I + K®X+ 

A(X-) = X-®K-'^ + I®X- (5.10) 

The operation A(-) is called the comultiplication. In the quantum group, the comultiplication does 
not commute with the exchange operator u defined hy aa ® h = h ® a. However, there is an operator 
R which satisfy the following 

Rl^{x)=a/\{x)R, for x e Ug^sh) (5.11) 
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Thus, the tensor product Vi V2 can be related to V2 ® Vi through the R matrix of the quantum 
group. 

In the Uq{sl2), the R matrix can be constructed in the operator formahsm. If operators X"^ and 
H satisfy the defining relations of the Uq{sl2), then the operator TZ defined by the following satisfy 
the intertwining relation (5.11) 

where exp^ x denotes the infinite series 

00 

exp,. = $:5-"^"-^^/^^ (5.13) 

n=0 \-"'^- 

The operator TZ is called the universal R matrix. 

The representation (5.8) of the braid group corresponds to the representation of the universal R 
matrix on the tensor product of two fundamental representations. 

We note that the universal R matrix can be constructed canonically through Drinfeld's quantum 
double construction (for instance, see [24]). This is similar to the Sugawara construction which derives 
the energy momentum tensor from the current operator. 
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A Commuting transfer matrices and the Yang-Beixter equations 



We show that if given two sets of the Boltzmann weights of the six vertex model satisfy the Yang- 
Baxter relation, then their transfer matrices commute. We consider three sets of the Boltzmann 
weights: {wi,W2,W3) = {a,b,c), {a ,b ,c), and (o ,b ,c ). Let us denote by r and r the transfer 
matrices constructed from the sets of the Boltzmann weights (a , 6 , c ) and (a ,b ,c ) , respectively. 
Then, we can show that if the three sets of the Boltzmann weights satisfy the Yang-Baxter equations 
given in §2.3, then the transfer matrices r and r commute. 

Let us now explicitly discuss the commutation relation. We first introduce the monodromy matrix. 
It is an N ranked tensor, whose {a, [3] elements are defined as follows 

{Ta,i3)l\]Z'^^ = ^ w{(3,hi\auC2)w{c2M\a2,c^)- ■ ■w{cN,hN\aN,Oi) (A.l) 

C2,...,Cjv 

The transfer matrix is given by the trace of the monodromy matrix 

r = tr(r)= Ta,a (A.2) 

a=l,2 

In terms of the matrix elements, we have 

a=l,2 



ai 



a2 



UN 



P. 





C2 


CN 






bi 


b2 


bN 



a 



Figure A.l: The matrix element (a, (3) of the monodromy matrix (Tq^^)^^''"'^^ 

Let us denote by T' and T" the monodromy matrices for the sets of the Boltzmann weights 
{a ,b ,c) and (a ,b ,c ), respectively. Wc consider the product of the matrix elements of the two 
monodromy matrices: p'^y S- ^^^^ entry of (ai, . . . , oat) and . . . , bN) of the product ^T^ ^ is 
given by 

\ ei,...,ejv 



(t' 

V 7l,7JV+l- 



|5l,<5A^+l ) 



6i,...,6jv 



6i,...,&jv 



ei,...,ejv 

= X] E w'{'ji,ei\ai,C2)w'{c2,e2\a2,C3)---w'{cN,eN\aN,lN+i) 

ei,...,ejv C2,...,cjxT 

X ^ w"{Si,bi\ei,d2)w"{d2,b2\e2,d3)---w"{dN,bN\eN,SN+i) 

6,2, ■■■,dN 

= E E E (^'(7i,ei|ai,C2)w"((5i,6i|ei,(i2)) • 

C2,---,cjv d2,-.;dN ei,...,eiv 

(w' {c2,e2\a2,C3)w" {d2,b2\e2,d3)^ ■ ■ ■ (w' {cN,eN\aN,lN+i)'w" {dN,bN\eN,SN+i)) 
= E S{aub^)]lX■S{a2MZ2■■■S{aNM)TS.T^ (^-4) 



C2,...,CN d2,...,dN 



Here S{aj,bj)'^J.''^J^^^ has been defined by 



Siaj, bjfj.'2^^[ = J2 w' {cj,ej\aj,Cj+i)w" idj,bj\ej,dj+i) 



(A.5) 
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We define the matrix element as follows 

Ml2=^idoMc<^,c,) (A.6) 

Here we assume that denotes the matrix element for column (co,(io) and row (ci,(ii) of the 

matrix M. Multiplying the matrix M to the product (i^'-y 5 ^^"^ ^PPtyiiig the Yang-Baxter relation 
N times, we can derive the following 

-'"i5o,<ii cijTjv+i-^dijdiV+i cijCjv+i-'dijdjv+i djVi^AT+i ' ) 



OAT 



ai 



02 



AT+l 



70 








J 7iv+i 70^ 


Cl 


C2 J 




,7Af+l 








! « 








) ( 






di 


6l 


62 


Sat 


di 

hi 


d2 1 


62 





70 1 Cl 



CAT 7A''+1 




Figure A. 2: Pictorial proof of the commutation relation: MT't" = T"t' M. Open and closed circles 
denote the Boltzmann weights w and w , respectively. The summation over variables ci, . . . , cjv and 

di, . . . , d^ is assumed. 

Let us briefly discuss the derivation of relation (A. 7). It is depicted in Fig. A. 2. In the first 
equality of Fig. 2.6, we have applied the Yang-Baxter relation formulated as follows 



(A.8) 



Here the symbol S' {aj,bj)'^J^'^^'^^^ has been defined by 



S' (oj , = E ^" (^J ' ' Cj+i)^' {dj , hj\ej , dj+i) 



(A.9) 



We also note that the LHS of (A. 7) corresponds to the sum 

E E M]^:alS{aiMr,\% ■ S{a„ h,)'^^, • • • 5(a^, h^)Z%X\ 

Cl,C2,—,CN dl,d2,---4N 



'ci,C2"dl,d2 



(A.IO) 



Let us consider the inverse of the matrix M 

V / di,d2 ^ / di,d2 

Multiplying the inverse to both hand sides of (A. 7) we have 

mt't"m-^ = t"t' 

Noting tr{MT'T"M-^) = tr{T'T"), we obtain the commutation relation of the transfer matrices 



(A.ll) 



(A.12) 



/ // // / 

r r = r T 



(A.13) 
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